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We present a covariant and supersymmetric theory of relativistic hydrodynamics in four- 
dimensional Minkowski space. 



Relativistic fiuid mechanics |l|] is generally simpler to formulate than its non-relativistic counterpart. As it is also 
believed to provide a more accurate description of hydrodynamical phenomena, much work has been invested in its 
development [gJjĮ. Recently, an interesting extension of the theory to include non-abelian charges and currents has 
been proposed One of the important aspects of this formalism is that it includes vorticity consistently at the 
lagrangean level, by developing a non-abelian generalization of the Clebsch decomposition of the vector conjugate to 

■ the current; for a review with many references, see [gj. 

£NJ \ In a related development, Jackiw and Polychronakos || have presented a supersymmetric theory of fluid mechanics 

■ in (2+l)-dimensional space-time. This model is rather special, as it descents from a supermembrane theory in (3+1) 
£Li! dimensions [@-^| . It results in a supersymmetric generalization of the non-relativistic planar Chaplygin gas []l0| . An 

interesting result obtained in |6| is, that the vorticity in the theory is generated by the fermion fields, rather then by 
the bosonic component of the fluid. In spite of these advances, so far a relativistic and supersymmetric theory of fluid 
mechanics in (3+1) dimensions is lacking. The present work contributes to filling this gap. 

The main result of this work is a supersymmetric component action, which we present both in its lagrangean and 
hamiltonian form, for a pseudo-classical conserved current made up of bosonic as well as fermionic contributions. 
To be more precise, the current is constructed from two real scalar fields (C, N) and a chiral fermion field V+Q : 

^' V ^gJČ) ^N+ l -G\C)^ +1 ^\ d-V = Q. (1) 

Here G(C) is some function of the real scalar field C which, for the purpose of constructing models, needs no further 
specification. However, under particular conditions discussed below the current allows an elegant and straightforward 
interpretation in terms of fluid flow. 

The fields (C, TV, iĮj±) and representing the conserved current are described by the supersymmetric lagrangean: 
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The infinitesimal supersymmetry transformations leaving the action invariant, parametrized by anti-commuting spinor 
parameters e±, are: 



S ■ SC = e+V>+ + e_^_, SN = iG(C) (e_V- - e+V+) , 

51> + = (ęC + iY)e-, 5^ = (ęC-if)e + , (3) 

<SV M = 2il + a^ v d y ^ + - 2ie^a fll/ d' y ip- 1 
Under these transformations the variation of the lagrangean is a totai divergence: 

SC = 9" [e+B +ll + e_S_ M ) , (4) 
where the vector-spinor fields B±^ are given, modulo eąuations of motion , by 
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x Our conventions for chiral spinors are such, that 7sV>± = iip± an d ^±75 = i^±; charge conjugations acts as ip± = Cip±, 
where %Į)± = iip^o. It should also be noted, that the euclidean 74 = 170 is hermitean, hence 70 is anti-hermitean. 
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The commutator of two supersymmetry transformations (^Į) closes as usual, modulo field eąuations. A complete and 
manifest off-shell formulation can be obtained from a superspace treatment. The relevant superspace action is defined 
in terms of a real vector superfield V with lowest components {C,iį)±,V ll ) 1 and a pair of conjugate chiral superficlds 
(A + , A_) with lowest component z such that Re z = N. The superfield expression for the lagrangean is 
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Here F(V) is an arbitrary real function of the superfield V, with lowest scalar component F(V)\g ±= o — F (C). 
Reduction of the superspace expression (|^) in terms of components leads, after elimination of a number of auxiliary 
fields and some rescaling of the scalars, to the expression (||) with the identification G(C) = F" (C). 
The field eąuations deri ved from the lagrangean (JŠJ) are: 



d-V = 0, 
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for the bosonic fields, and 



G(c) ^± = --g'(c) (ęc t if) V>± 
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(8) 



for the fermionic ones. The first two eąuations (0) indeed reproduce eąuation ([!]). 

Because of translation invariance and supersymmetry, the theory described by C conserves four-momentum and 
supercharge. The corresponding currents are provided by the energy-momentum tensor and the supercurrents. For 
the energy-momentum tensor we use the symmetrized version 
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It is obtained from the non-symmetric translational Noether currents O^j, by the addition of an improvement term: 
Tįi, v = + tipu, where 
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The field eąuations (ĮŽĮ) and (||) imply <9 M T M „ = 0, from which the conservation of four-momentum = J d 3 xT fl0 
follows. The supercurrents are obtained directly by Noether 's procedure and read (in an on-shell version) 
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As on shell d ■ S± = 0, the conservation of the supercharges Q± = J d 3 xS±o follows. 

The manifestly covariant lagrangean description of the theory has an eąuivalent canonical formulation in terms of 
a hamiltonian with a corresponding bracket structure. First, the canonical momentą are defined by 
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The inclusion of the 70 in the definition of the fermionic momentą is motivated by the general reąuirement that the 
charge conjugation properties of the momentą reflect those of the fermion variables themselves. Furthermore note, 
that Vq is a canonical momentum, whereas the 3-vector field V is an auxiliary field, which can be eliminated by its 
algebraic field eąuation; in particular in the following we use the identifications 



Va = n N j 



V = 
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As usual, the fermionic momentą are not independend, and the theory possesses the second-class constraints: 



(13) 
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As a result, the dynamics of the theory is generated by a hamiltonian — obtained from £ by a straightforward Legendre 
transformation — and a set of Dirac-Poisson brackets derived by reduction of the standard Poisson brackets to the 
physical shell (ĮliĮ) in the full phase-space. In summary, one first constructs a hamiltonian density 
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where it is to be noted, that f = 7 ■ V now denotes a 3-dimensional contraction. Subseąuently the field eąuations 
are obtained from the hamiltonian H = J d 3 xH by the Dirac-Poisson brackets 



d A = {A,H}% 
defined in terms of the elementary brackets/anti-brackets 

{C(r, t), n c (r>, t)}* = S 3 (r - r'), {N (r, t), n N (r>, t)}* = S 3 (r - r'), 

{n c (r, t), $ ± (r', t)}* = ^±(r, t) S 3 (r - r'), {V±(r, t), n c (r', t)}* = -^7* V±(r, i) <5 3 (r - r'), 



(16) 



2G(C) 

ftMr.t),^,*)}' = 7 ° 2 ( įĮg 5) ^(r-rO- 
In particular we note, that the eąuation 

^JV = {tTJV, H e ff}* = V • 
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after the identification (^ŠĮ) is eąuivalent with Vq = V -V , ov d-V — 0. A somewhat tedious calculation shows that (ĮT^) 
and ( ĮT^ ) indeed reproduce all covariant field eąuations (0) and (|Į). We can now construct the canonical expressions 
for the conserved four-momentum and the supercharges; for the energy-momentum vector we find the results 
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The canonical expressions for the supercharges are 
Q+= f d 3 r 
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Like the hamiltonian generates the time-evolution, the supercharges generate the supersymmetry transformations: 
the results (|Į) are reproduced by the brackets 
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The supercharges satisfy the standard super-Poincarė algebra, in particular {<9±,<3t}* = ^P- Note that the su- 
persymmetry transformation rule implies o;„ = {V fM ,Q±}* — 2i a fJiI/ d 1, '>jj±, where uj^ is a conserved fermionic current: 
d-LJ = 0, as expected from supersymmetry and the conservation law ol V^. It is also of interest to discuss the brackets 
of the current components among themselves. By applying the identifications (Įl^) one finds the non-trivial results 

{F (r ; i),f(r',t)}*^^V r( 5 3 (r-r'), and {^(r,t), V^r'^)}* = j +1 °a^- č 3 (r~r'). (22) 

Note that the charge-conjugation properties of the spinors guarantee that ip^^aijip- = iip_<JijiĮ)- = iip^_aįjip + , as 
in our conventions <r%j is anti-hermitean. 

For easy comparison we have performed the canonical analysis in terms of the ficlds related by the linear supersym- 
metry transformations (|Į), at the price of dealing with off-diagonal terms in the Dirac-Poisson brackets (0). Observe, 
however, that the brackets can be diagonalized by a field redefmition 

*± = y/G{Č)il>±. (23) 
Indeed, in terms of the new fermion fields the brackets read 

{C(r, t),7T C (r', t)}* = <5 3 (r - r'), {N (t, t), ir N (r>, t)}* = <5 3 (r - r'), 



{7r c (r,i),# ± (r',č)}* ={$ ± (r,i),7T C (r',t)}* =0, (24) 

|$ ± (r,t),$ T (r',t)}* =i/(l T75 ),5 3 (r-r'), 

but of course the supersymmetry transformations (||) now become non-linear. We do not present here the new 
expressions for the energy-momentum and supercurrents; however, the conserved current ([į]) is of particular interest 
and we observe that after the field redefmition ( p3| ) its space components read 

V = —L- V7V+ lG '^ l» + 7$ . (25) 
G(C) 2{G{C)) 2 +1 ~ K ' 

To complete its hydrodynamical interpretation, we relate the fields in our model to the fluid density p and velocity 
u p. First we consider the bosonic rcduction, obtained by reąuiring the fermion field to vanish: $± = 0. The field 
eąuations (ĮTĮ) then reduce to 
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V^pu^-^—d^N, vl = - P \ d-v = o, 
G{C)uC= l -G'{C)(-(dCf + p 2 ). 



(26) 



In the present context the current-conservation condition is to be interpreted as the eąuation of continuity in hydro- 



dynamics, whilst the last eąuation (26) can be used to express the density p in terms of the field C. The first eąuation 
( p6| ) then implies that the velocity field for fixed p (i. e. for fbced C) depends only on one degree of freedom, N. This is 
also the situtation encountered in potential flow, when u M = <9 M 6>, (dO) 2 = — 1. Therefore our bosonic model describes 
potential flow with vanishing vorticity, providcd 
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Such a relation is consistently realized if either d^p ~ 8^9 ~ V lll or else if N can be expressed entirely in terms of 
9, independend of p: dN/dp = 0. In the first case d^C ~ Vį t , and the bosonic part of the energy-momentum tensor 
manifestly takes the perfect fluid form; indeed, if d^C — nu^, then 

Tffi = P9^ + (p + £)«„«*.. (28) 



4 



where the pressure and energy density are given byĘ p = e = | G(C) (k 2 + p 2 ). The case dN/dp — O is realized 
if JV = 2m#, p = 2m/G(C), with m a constant of proportionality It follows that N must satisfy the condition 
(dN) 2 = —Am 2 = constant, whilst combination of these conditions with ( f26| ) gives 

^|^^(v^(č)^)=2m 2 . (29) 

For example, if G(C) = l/C, we find the imaginary-mass Klein-Gordon eąuation 

(□ + m 2 ) VČ = 0, (30) 

which allows non-trivial stationary solutions for the density p = 2mC. The more general Ansatz G (C) — C p leads to 
the non-linear relation 

□ C^' 2 = m 2 p(p + 2) C- ( - 3p + 2 ^ 2 . (31) 
If the fermion fields are switched on, the potential character of the fiow disappears; indeed, as 

v, =p-,= o{č) W + ^li? * +7 ^' (32) 



the velocity is now parametrized in terms of one bosonic and two fermionic degrees of freedom. As eq.(|22j) shows, the 
vorticity is non-zero as a result of the fermionic contribution indeed. In terms of the redefincd fields <Į> these bracket 
relations take the form 

{V-(r,i),V,(r',t)}* = (|§į§jj) 2 *+7%<i>- č 3 (r-r'). (33) 

In particular, for the case G(C) — l/C the coefficients of the bi-fermion term on the r. h. s. of eqs.(|32]) and ( ĮŠŠĮ ) reduce 
to constants, and the current is a linear combination of free boson and fermion currents. 

The models constructed here can be generalized to include vorticity from bosonic potentials, by a supersymmetric 
extension of the Clebsch decomposition of V^; this is explained elsewhere jll]]. Consistent extensions of the models at 
the ąuantum level are of potential interest in cosmology, where they could provide an effective description of an early 
supersymmetric phase of the universe, and in condensed matter physics, where they might apply to ąuantum fiuids 
like 3 He- 4 He mbctures, up to effects proportional to the mass-differences of these isotopes. 
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